Coexistence of dimerization and long-range magnetic order in the frustrated spin 
chain system LiCu202: inelastic light scattering study 
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Raman scattering studies of the frustrated spin chain system LiCu202 are reported. Two transitions 
into a magnetically ordered phase (taken place at temperatures ^ 9 K and ^ 24 K) have been 
confirmed from the analysis of optical properties of the samples. Interestingly, two different magnetic 
excitations, seen at 100 and 110 cm~^ in the magnetically ordered phase superimpose each other 
independently, indicating a coherent coexistence of long-range magnetic order and dimerization. 
The observed phenomenon is attributed to magnetostructural peculiarities of LiCu2 02 leading to 
the intrinsic presence of nonmagnetic impurities on a nanometer scale. Furthermore, magnetic 
impurities play a significant role in driving the transition from an incommensurate state to a Neel 
ordered one at 9 K. 



Theoretical and experimental investigations of low- 
dimensional spin systems with frustration and dimeriza- 
tion have been boosted by the discovery of the inorganic 
spin-Peierls system CuGeOa. [1] Such a frustrated spin 
chain system shows a rich phase diagram, including both 
gapped and gapless phases [2] , as well as disorder- induced 
long-range ordering, coexisting with a dimerized state. 
[3] Variational calculations suggest that the resonating- 
valence-bond character of spin correlations at short dis- 
tances can be responsible for the enhancement of anti- 
ferromagnetic (AFM) correlations near vacancies [4,5], 
that eventually results in long-range AFM order, locally 
coexisting with the disordered dimerized phase. Such 
a coexistence was observed in some other doped systems 
with reduced dimensionality (for instance, in highly hole- 
doped chain material Sro.73Cu02 [6]) and appears to be a 
fundamental property of low-dimensional quantum spin 
systems. 

LiCu202 can be regarded as a reahzation of a S—1/2 
spin chain with competing nearest- and next-nearest- 
neighbor interactions. [7-11] This compound has an or- 
thorhombic crystal structure of a space group Pnma with 
the lattice parameters a=5.72 A, 6=2.86 A and c=12.4 A. 
[7] There are monovalent and bivalent copper ions in the 
unit cell. Magnetic Cu^+ ions form a double-chain along 
the b axis which is separated from each other by both 
Li ions and planes with nonmagnetic Cu"*" ions. The 
two separate Cu chains within the double-chain struc- 
ture are coupled via a 90° oxygen bond along the c axis. 
At elevated temperatures high-field electron spin reso- 
nance (ESR) gives evidence for a spin singlet state with 
spin gap of A ^ 72 K. [10] Interestingly, upon coohng a 
spin singlet state transits into a long-range ordered state 
with helimagnetic structure at ~ 24 K. [11] Some 
bulk measurements point to the presence of a second low 



temperature transition with a coUinear AFM structure 
at ~ 9 K. [10] Both transitions are attributed to 
an intrinsic non-stoichiometry and the effect of nonmag- 
netic and/or magnetic impurities. [10,11] However, the 
exact origin is not yet clear. In addition, with decreas- 
ing temperature a reduction of orthorhombic strain has 
been observed. [8] These magnetostructural peculiarities 
of LiCu202 provide a good opportunity to investigate the 
influence of a variety of magnetic and nonmagnetic im- 
purities on magnetic property of a frustrated spin chain 
system. 

Raman spectroscopy has proven to be an extremely 
powerful technique to probe magnetic excitations and 
spin-lattice interactions in a low-dimensional spin sys- 
tem with unprecedented precision. [2] The main moti- 
vation of this work is using Raman spectroscopy tech- 
nique to study the nature of ground state and the mag- 
netostructural peculiarities in quantum AFM LiCu202. 
In the following we will show the coherent coexistence of 
long-range ordered states with dimerized state induced 
by nano-scale defects. Further, we propose a possible 
mechanism for a disorder-induced long-range ordering in 
the studied system. 

The single crystals were grown using a self-flux method 
and characterized by a microstructural analysis and ther- 
modynamic measurements as described in Ref. [10]. The 
single crystals are microscopically twinned and contain 
LiCuO-impurity phase. Raman spectra were measured 
in a quasi-backscattering geometry with the excitation 
line A = 514.5 of an Ar+ laser with the power of 10 mW 
and were analyzed by a DILOR-XY spectrometer and a 
nitrogen cooled charge-coupled device detector. 

Figure 1 displays Raman spectra in parallel (xx) and 
crossed (xy) polarizations at 3 K as well as at 5 K and 
room temperature in xx polarization. Here the x axis is 
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an arbitrary direction in the ab plane. Raman spectra in 
xy polarization exhibit the same behavior as those in xx 
polarization with weaker intensity due to a twinning of 
the single crystal. [10] Thus, the observed spectra can be 
regarded as an average of all ab plane polarizations. Sub- 
tracting the acoustic (Bi„ + B2„ + B3„) modes the factor 
group analysis of the space group Pnma yields the follow- 
ing Raman- and infrared-active modes; 10 Ag(aa,bb,cc) 
+ 5 Big(ab) + 10 B2<,(ac) + 5 B3g(bc) 5 A„ + 9 Bi„ 
-I- 4 B2„ -|- 9 B3„. At room temperature we observe 12 
Raman-active modes out of the maximally expected 15 
modes in ab plane polarizations. A broad band extend- 
ing from 850 to 1300 cm~^ is an overtone feature of the 
first order signals between 400 and 650 cm~^. This might 
be due to strong anharmonic lattice interactions and/or 
resonant scattering. 
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FIG. 1. Raman spectra of LiCu202 in parallel (xx) and 
crossed (xy) polarizations at 3 K (upper panel) as well as 
at low and room temperature in parallel polarization (lower 
panel). The arrows indicate additional phonon and magnetic 
signals which appear at low temperature. Inset: Temperature 
dependence of the 494- and 562- cm~^ phonon modes. 

With decreasing temperature both new phonon modes 
and magnetic continua show up in an intriguing man- 
ner. First, we focus on phonon anomalies which reflect 
changes of the local symmetry and bond strengths. Upon 
cooling new phonon modes at 475 cm~^ and 518 cm^^ 
appear below 55 K and 200 K, respectively. The former 
temperature corresponds to roughly the exchange con- 
stant of J ~ 66 K [11] while the latter one is related 
to the appearance of the spin gap. [10] Note here that 
the number of the observed 14 modes in the low tem- 
perature region is still smaller than the predicted one 
(15 modes), suggestive of no reduction of crystal symme- 
try at low temperatures. In addition, almost all phonon 
modes become sharp and intense while the Cu-0 bending 



and stretching modes show significant changes as the in- 
set of Fig. 1 displays. Upon cooling the 494-cm~^ and the 
562-cm~^ modes undergo an appreciable hardening of 5-8 
cm~^ and then saturate around T = 55 K. Upon further 
cooling they soften slightly while showing a tiny jump 
at 50 K. Noticeably, a softening takes place around tem- 
perature of the exchange constant J 66 K. [11] Thus, 
a softening below 50 K can be attributed to a renor- 
malization of phonon energy via spin-phonon coupling, 
similar to the alternating spin chain system (VO)2P207. 
[12,13] The smallness of the observed softening is related 
to the big difference in energy scale between J ~ 66 K 
and the optical phonon energy of 700 K which me- 
diates the exchange paths along the double-chain. High- 
resolution x-ray diffraction measurements [8] show a sub- 
stantial increase of orthorhombic strain [a — b)/(a + b) 
upon heating. Transition metal oxides may show a de- 
crease of the orthorohombicity with increasing tempera- 
ture, since thermally activated lattice vibrations reduce 
the strain. Normally, phonon frequencies shift to higher 
energy, linewidths broaden, and integrated intensity of 
some phonon modes decreases. The studied system fol- 
lows the behavior expected for usual systems despite the 
opposite behavior of the orthorhombicity as a function of 
temperature. To clarify the seemingly inconsistency be- 
tween local and bulk properties detailed studies of struc- 
ture as a function of temperature are needed. 

We will turn now to distinctive features of magnetic 
excitations observed at low temperatures and low fre- 
quencies. As Fig. 2 displays, at 3 K a broad asymmetric 
continuum around 100 cm~^ is seen. With increasing 
temperature the continuum shifts to lower energy while 
damping increases. At the same time, a quasielastic scat- 
tering response and a weak continuum around 110 cm~^ 
develop. Although the crystals are twinned, the dynam- 
ics of low-dimensional spin system can be well resolved by 
Raman spectroscopy. This is because magnetic Raman 
scattering of low-dimensional systems contributes selec- 
tively to the polarization direction in which the incident 
and scattered light is parallel to the dominant exchange 
paths. Thus, twinning does not add an essential diffi- 
culty to interpreting magnetic signals. To clarify the rela- 
tion between the magnetic excitations and the structure 
of magnetic ordering observed in the bulk material we 
will distinguish three phases following Ref. [10]; I: mag- 
netic ordered phase (T < 9 K, presumably with coUinear 
AFM structure), II: helimagnetic ordered phase [14] (9 
K < T < 23 K), and III: dimerized phase (T > 23 K). 

In phase I the broad continuum extending from 40 
cm^i to 130 cm ^ is observed. This continuum per- 
sists as a wing feature of quasielastic scattering up to 
T=21 K (2.3 Tea). This feature is typical for two-magnon 
(2M) scattering originating from double spin-flip pro- 
cesses via the exchange mechanism in antiferromagnets 
with coUinear structure. [15] Therefore, we identify phase 
I to be Neel ordered with T^^^T^^ = 9 K. 



2 



The evolution of the 2M spectrum reflects mainly the 
temperature dependence of short-wavelength magnon en- 
ergies and lifetimes. [16] Thus, the persistence of 2M scat- 
tering to several Tjv can be interpreted in terms of the 
presence of short-range magnetic fluctuations damped 
by thermal fluctuations. In Fig. 3 its temperature 
dependence of the normalized 2M frequency and the 
full width at half-maximum is presented together with 
higher dimensional results. [16,17] Magnon-pair energies 
of LiCu202 (ID S=l/2) are renormalized only by 3% 
at Tjv. In contrast, the magnon-pair spectral weight is 
renormalized by 25 % at Tn for 3D S=l/2 systems and 
by less than 5% for 2D S=l systems. [16,17] The damping 
does scarcely take place at Tm for LiCu202. However, it 
strongly increases as the dimension and spin number in- 
crease as the right panel of Fig. 3 displays. The higher di- 
mensionality and spin number are, the larger are changes 
of spectral weights at an energy scale comparable to the 
Neel temperature. This is related to the fact that the 
Neel temperature is not an appropriate energy scale for 
magnetic excitations in low-dimensional systems. Com- 
pared to higher dimensional systems, the robustness of 
spin-fluctuation dynamics at the energy scale of Tjv con- 
firms the low-dimensional character of the studied sys- 
tem. 
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FIG. 2. Low frequency Raman scattering of LiCu202. 
Quasielastic scattering and two different kinds of two-magnon 
continua have been observed. 

The frequency of the 2M peak allows an estimate of the 
unrenormalized AFM exchange constant between copper 
spins. The exact determination of the exchange constant 
in the studied system is impossible because no theoreti- 
cal calculation for a double-chain system is known. How- 
ever, one can make a reasonable estimation using not 
too stringent assumptions. Because of frustrations, the 
magnetic behavior of LiCu2 02 lies between ID and 2D 
antiferromagnet. In the 2D case, the peak energy of 2.7 J 
imposes a lower boundary of the exchange constant. [18] 



In frustrated chain systems the one-magnon dispersion is 
given by uj{k) = Jil — acos(kd) /2) with the frustration 
parameter a and the chain repeat vector d. [19] In the 
noninteracting case 2M scattering is given by twice the 
magnon density of states. Thus, the peak energy corre- 
sponds to J [2 + a). If we assume a renormalization of the 
peak energy due to magnon-magnon interactions by the 
order of J a, then an upper boundary of the exchange 
constant is roughly given by 2J. To conclude, one ob- 
tains 52 K < J < 70 K from the peak energy of 140 K. 
This value encompasses J = 66 K obtained by a fit of a 
frustrated chain model to the static susceptibihty. [11] 
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FIG. 3. Comparison of renormalized frequency (left panel) 
and damping (right panel) of two-magnon continuum as a 
function of spin number and dimensionality. The higher di- 
mensional data are taken from Ref. 16. The dashed vertical 
line marks a Neel temperature. 



In the following we will discuss the phase II and HI. 
Quasi-elastic scattering develops in the respective tem- 
perature interval. In addition, a weak continuum ex- 
tending from 102 cm^^ to 120 cm^^ becomes clearly dis- 
tinguishable from the broad 2M continuum. This con- 
tinuum persists up to 200 K, that is, into the dimerized 
phase HI. We thus assign the observed signal to a 2M 
continuum which corresponds to a double spin-fiip pro- 
cess of two singlets into a higher singlet state. [2] Since 
the onset of the continuum corresponds to twice the spin 
gap, one obtains a spin gap of A w 73 K. This value is in 
excellent agreement with A « 72 K from high-field ESR 
measurements. [10] 

However, we find no evidence for the presence of ad- 
ditional magnetic excitations in the magnetic long-range 
ordered state II. The silence of magnetic Raman scat- 
tering seems to be related to the incommensurate heli- 
magnetic structure suggested in Ref. [11]. In contrast 
to the 2M continuum from the Neel ordered state, the 
2M continuum from singlet states shows no noticeable 
change in peak position as well as intensity as a func- 
tion of temperature. This is due to the narrow band- 
width of triplet excitations compared to the magnitude 
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of the spin gap. Moreover, two 2M continua are simply 
superposed to each other independently. This indicates 
that in phase II short-range singlet correlations and the 
helimagnetic ordered state coexist coherently. Further- 
more, the simultaneous observation of magnetic excita- 
tions from Neel order and singlet correlations in phase I 
can be interpreted as the coherent coexistence of an an- 
tiferromagnetic long-range ordered and a dimerized state 
which are mutually exclusive. Note here that the dimer- 
ized state is robust in the course of the evolution of mag- 
netic structures and their corresponding magnetic cor- 
relations. As mentioned above, such a coexistence has 
been reported in the impurity-doped spin-Peierls system 
CuGeOs [3] and the highly hole-doped quasi- ID cuprate 
Sro.73Cu02. [6] The introduction of nonmagnetic impu- 
rities to dimerized states produces a local moment by 
breaking up dimers. With aid of higher dimensional in- 
teractions long-range ordering occurs at low temperature. 
[4] This disorder-order transition scenario might be ap- 
plicable to the studied system as the dimerized phase of 
LiCu202 is intrinsically contaminated by a LiCuO impu- 
rity phase which has copper ions in a nonmagnetic Cu"*" 
oxidation state. The LiCuO phase is estimated to be 
less than 10% of the total volume. It is arranged almost 
regularly in the form of platelets with the dimension of 
100 X 7 X 100 nm^ (Ref. [10]). It is thus natural to expect 
that similar to a single-site dopant, such nanostructural 
nonmagnetic inclusions break Cu^+-Cu^^ bonding along 
the chains, enhance three-dimensional interactions and 
eventually promote AFM long-range order at low tem- 
peratures. 

However, within this mechanism it is difficult to cap- 
ture the following features. First, the presence of a he- 
limagnetic state is not explainable because for the in- 
troduction of dopants the interaction between local mo- 
ments has an alternating sign. That is, its sign relies on 
whether the two local moments belong to the same or to 
opposite sublattices. As a result, frustrations are lifted 
and an AFM ordering is favored. Second, the transition 
of an incommensurate phase into the Neel ordered one at 
9 K is driven by the presence of the magnetic Li2Cu02- 
impurity. [11] Third, the 2M signal from a long-range or- 
dered state is much stronger than that from a dimerized 
phase (see Fig. 2). All these observations indicate the 
predominance of a long-range ordering over a dimeriza- 
tion as well as the dependence of the magnetic structure 
on the kind of impurities, signalling the significant role 
of magnetic impurities. In the following we present a 
possible scenario. 

Upon cooling dilute magnetic impurities which lie most 
probably between spin chains undergo long-range order- 
ing via 3D interactions. Then, the long-range ordered 
impurities can be regarded as effective magnetic fields in- 
ducing long-range ordering in spin chains. The observed 
helimagnetic state is due to the compromise of AFM or- 
der with strong in-chain frustration. Furthermore, two 



successive transitions at 22.5 K and 24.2 K [10] reflect 
an intriguing interplay of magnetic impurities between 
chains and nonmagnetic in-chain impurities on the back- 
ground of the incommensurate spin structure. Upon fur- 
ther cooling, depending on the presence of the Li2Cu02- 
impurity, the Neel ordered phase will appear at about 9 
K. This indicates that the coUinear magnetic structure 
of the Li2Cu02-impurity promotes the transition of the 
incommensurate phase to the commensurate one. 
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FIG. 4. (a) Comparison of the intensity of quasielastic scat- 
tering (full rectangle) and ESR (open triangle) to static sus- 
ceptibility (solid line) taken from Ref. 10. (b) Mapping of 
the magnetic specific heat derived from quasielastic scatter- 
ing (full rectangle) on the magnetic specific heat obtained by 
subtracting the phonon contribution (solid line). 

Finally, we will discuss the quasielastic Raman re- 
sponse observed above 10 K. Our experimental setup 
is carefully adjusted to suppress Rayleigh scattering. 
Therefore, the observed intensity is intrinsic and origins 
from fluctuations of the energy density of the system. 
Its presence is a remarkable feature of low-dimensional 
system with strong spin-phonon coupling. [2] 

According to the theory of Reiter [20] and Halley [21], 
the scattering intensity is given by the Fourier compo- 
nent of the correlation function of the magnetic energy 
density: I{uj) oc J^^dte^"^* < E{k,t)E*{k,t) >, where 
E{k,t) is the magnetic energy density. In a hydrody- 
namic assumption for the correlation function in the high 
temperature limit [22] the above equation is simplified to 
the Lorentzian profile oc ^r^^jy^^yi^ where k is the 
scattering wave vector, Dt the thermal diffusion con- 
stant, and Cm the magnetic specific heat. In this case, 
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magnetic specific heat is proportional to the integrated 
intensity divided by T^. Using this relation one can 
map scattering intensity on thermodynamic quantities. 
Figure 4(a) displays integrated intensity of quasiclastic 
scattering together with ESR intensity and static sus- 
ceptibility as a function of temperature. There exists an 
overall good correspondence between integrated ESR and 
quasielastic intensity. However, they show a substantial 
deviation from the static susceptibility below 45 K. Here 
note that ESR probes excitations within excited states 
which originate from the low-dimensional character of in- 
teractions in LiCu202 as the drastic drop of its intensity 
below 45 K shows. Thus, one can see that ESR and Ra- 
man spectroscopy are selectively sensitive to disordered 
short-range correlations. In contrast, the static suscepti- 
bility becomes dominated by 3D short-range-order anti- 
ferromagnetic correlations below 45 K. In Fig. 4(b) the 
magnetic specific heat derived from the quasielastic scat- 
tering is shown together with the magnetic part of spe- 
cific heat obtained after subtracting a calculated phonon 
contribution from the measured specific heat. In the tem- 
perature interval 12K<T<40K there is a reasonable 
matching between them. The difference above 40 K can 
be attributed to an overestimation of the phonon con- 
tribution to the specific heat at high temperatures by 
choosing only one Debye function. The discrepancy be- 
low 12 K comes from the suppression of fluctuations of 
spin energy density in the Neel ordered phase. 

A maximum of C'm can also provide an information on 
the exchange constant. In the case of ID AF chain the 
magnetic specific heat has a broad maximum at /c^T « 
0.481J. The frustration shifts the maximum of Cm to 
lower temperature, fc^T ~ 0.38J. [23] This results in 
J ~ 60 K which is consistent with the value obtained 
from the peak position of the 2M scattering as well as 
from the static susceptibility. 

In summary, a Raman scattering study of the frus- 
trated spin chain compound LiCu202 unveils a coexis- 
tence of dimerized correlations and long-range order as 
a coherent superimposition of a dimerized state with a 
spin gap A ~ 73 K on a two-magnon continuum from a 
Neel ordered state shows. This is attributed to nano-scale 
sized nonmagnetic defects. Moreover, the transition of a 
helimagnetic state into a Nccl ordered one suggests a sig- 
nificant role of magnetic impurities in an incommensurate 
magnetic structure as a driving impetus to a commensu- 
rate structure. We hope that our experimental observa- 
tions will stimulate further theoretical investigations de- 
voted to understanding magnetic properties of frustrated 
spin chain systems including magnetic/nonmangetic de- 
fects. 

The work was supported by DFG/SPP 1073, IN- 
TAS 01-278 and NATO Collaborative Linkage Grant 
PST.CLG.977766. 



[1] M. Hase et al, Phys. Rev. Lett. 70, 3651 (1993). 
[2] P. Lemmens et al., Phys. Rep. 375, 1 (2003) and refer- 
ence's therein. 

[3] L. P. Regnault et al., Europhys. Lett. 32, 579 (1995). 

[4] G. B. Martins et al, Phys. Rev. Lett. 78, 3563 (1997). 

[5] E. S. S0rensen et al., Phys. Rev. B. 58, R14701 (1998). 

[6] G. I. Meijer et al, Phys. Rev. B 60, 9260 (1999). 

[7] R. Bergcr et al., J. Alloys Compd. 184, 315 (1992). 

[8] B. Roessh et al, Physica B 296, 306 (2001). 

[9] Vorotynov et al., JETP 86, 1020 (1998). 
[10] S. Zvyagin et al, Phys. Rev. B 66, 064424 (2002). 
[11] T. Masuda et al, cond-niat/0310126 (2003). 
[12] M. Grove et al, Phys. Rev. B 61, 6126 (2000). 
[13] U. Kuhlmann , et al., Phys. Rev. B 66, 064420 (2002). 
[14] Strictly speaking, the phase transition into a heUmag- 
netic phase around 23 K consists of two successive tran- 
sitions at 22.5 K and 24.2 K [10]. The possible origin will 
be discussed below. 
[15] P. A. Fleury et al, Phys. Rev. 166, 166 (1968). 
[16] M.G. Cottam and D.J. Lockwood, Light Scattering in 

Magnetic Solids, John Wiley & and Sons (1986). 
[17] P. A. Fleury et al, Phys. Rev. Lett. 24, 1346 (1970). 
[18] K. B. Lyons et al, Phys. Rev. B 37, 2353 (1988). 
[19] D. A. Tennant et al., Phys. Rev. B 67, 054414 (2003). 
[20] G. F. Reiter, Phys. Rev. B 13, 169 (1976). 
[21] J. W. Halley, Phys. Rev. Lett. 41, 1605 (1978). 
[22] B. I. Halperin et al, Phys. Rev. 188, 898 (1969). 
[23] H. Kuroe et al., Phys. Rev. B 55, 409 (1997). 



5 



